arXiv:l502.01224vl [astro-ph.HE] 4 Feb 2015 


Astronomy & Astrophysics manuscript no. castelletti-astroph ©ESO 2015 

February 5, 2015 


A complete radio study of SNR G15.4+0.1 from new GMRT 

observations 

L. Supan 1 , G. Castelletti 1 , B. C. Joshi 2 , M. P. Surnis 2 , and D. Supanitsky 1 

1 Instituto de Astronomfa y Flsica del Espacio (IAFE, CONICET-UBA), CC 67, Sue. 28, 1428 Buenos Aires, Argentina 
e-mail: lsupan@iafe.uba.ar 

2 National Centre for Radio Astrophysics (NCRA-TIFR), Post Bag No 3, Ganeshkhind, Pune - 411007, India 
Received <date>; Accepted <date> 


ABSTRACT 

Aims. The supernova remnant (SNR) G15.4+0.1 is considered to be the possible counterpart of the y-ray source HESS J1818-154. 
With the goal of getting a complete view of this remnant and understanding the nature of the y-ray flux, we conducted a detailed 
radio study that includes the search for pulsations and a model of the broadband emission for the SNR G15.4+0.1/HESS J1818-154 
system. 

Methods. Low-frequency imaging at 624 MHz and pulsar observations at 624 and 1404 MHz towards G15.4+0.1 were carried out 
with the Giant Metrewave Radio Telescope (GMRT). We correlated the new radio data with observations of the source at X-ray and 
infrared wavelengths from XMM-Newton and Herschel observatories, respectively. To characterize the neutral hydrogen (HI) medium 
towards G15.4+0.1, we used data from the Southern Galactic Plane Survey. We modelled the spectral energy distribution (SED) using 
both hadronic and leptonic scenarios. 

Results. From the combination of the new GMRT observations with existing data, we derived a continuum spectral index tt=-0.62 ± 
0.03 for the whole remnant. The local synchrotron spectra of G15.4+0.1, calculated from the combination of the GMRT data with 
330 MHz observations from the Very Large Array, tends to be flatter in the central part of the remnant, accompanying the region 
where the blast wave is impinging molecular gas. No spectral index trace was found indicating the radio counterpart to the pulsar 
wind nebula proposed from X-ray observations. In addition, the search for radio pulsations yielded negative results. Emission at far- 
infrared wavelengths is observed in the region where the SNR shock is interacting with dense molecular clumps. We also identified HI 
features forming a shell that wraps most of the outer border of G15.4+0.1. Characteristic parameters were estimated for the shocked 
HI gas. We found that either a purely hadronic or leptonic model is compatible with the broadband emission known so far. 

Key words. ISM: individual objects: SNR G15.4+0.1, HESS J1818-154-ISM: supernova remnants-Gamma rays: ISM-Radio con¬ 
tinuum: ISM 


1. Introduction 

The radio source G15.4+0.1 was first identified as a super¬ 
nova remnant (SNR) during the 90 cm multiconfiguration Very 
LargeAmwJ i VLA}_surve^_o^he_Galactic_glan^JBrogan 
st al. 120060 . This remnant is found to positionally coincide 
with extended emission (~8.5 aremin in size) detected in the 
TeV y-ray band with the High Energy Stereoscopic System 
(H.E.S.S.) dHofverberg et al.ll201 ll) . The TeV source, catalogued 
as HESS J1818-154, is one of the faintest HESS sources de¬ 
tected to date. On the basis of the morphological correspondence 
between the brightest hotspot of HESS J1818—154 and the inner 
part of G15.4+0.1, it has been proposed that a pulsar wind nebula 
(PWN) powered by a still undetec ted rotating neutron star could 
be responsible for the TeV y rays dHofverberg et al.ll201 ll) . 

Recentlv JCastelletti et al.ld2013l) (C13) found morphological 
and kinematical evidence for shocked CO gas in a dense molec¬ 
ular cloud located in the foreground of the SNR G15.4+0.1. 
Because of the location of the cloud (partially overlapping the 
HESS source) and the signatures of interaction between the 
molecular gas and the remnant, the authors argued that the 
gamma radiation might be arising from hadronic interactions 
produced by the passage of the supernova shock front through 
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the molecular cloud. In that paper, based on 21 cm HI absorp¬ 
tion measurements along the line of sight to G15.4+0.1, those 
authors also redetermined the distance to this remnant placing it 
at4.8 ± l i 0k£c J jnsteadofthe_10 - ±^k£c_estimatedb^Hofver- 
berg et al. (12011 ) using the uncertain X-A relation for supernova 
remnants. 

Later on, lAbramowski et all 12014}) used the XMM-Newton 
telescope to investigate the X-ray counterpart to the y-ray source 
HESS J1818-154. The new X-ray image of the remnant shows 
extended emission in the interior of the radio SNR shell, in 
spatial coincidence with the discovered HESS source and the 
southern part of the molecular cloud. Although the X-ray emis¬ 
sion is fitted satisfactorily by both the thermal and non-thermal 
models, these authors interpreted the X-rays as emission from a 
PWN. In addition, the X-ray observations lead to the detection 
of five point-like sources in the region of the y-ray emission. 
From their X-ray properties, only two were thought to be galactic 
sources that could be candidates for the central compact object 
of the SNR. The analysis of periodicity in the signal from these 
sources, however, has not shown evidence of X-ray pulsations. 

In this paper, based on new data obtained with the Giant Me¬ 
trewave Radio Telescope (GMRT), we provide the hist compre¬ 
hensive analysis of the morphological and spectral radio proper¬ 
ties of SNR G15.4+0.1, which includes the results of the search 
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for radio pulsations in the region of the source HESS J1818—154. 
In addition, we report on the study of the spatial correspondences 
between the radio, infrared, and X-ray emission bands, and the 
analysis of the neutral hydrogen distribution in the environment 
of the remnant. On the basis of the flux measurements from radio 
to y ray, we discuss scenarios in which the TeV emission from 
this source originates from either hadrons interacting with dense 
interstellar material or leptonic emission. 

2. Observations and imaging 

2.1. Radio continuum and neutral hydrogen observations 

The SNR G15.4+0.1 was observed using the GMRT at 624 MHz 
on June 17 and 18, 2012 (Project code 22-015). The project also 
includes GMRT observations of the remnant at 1420 MHz. The 
image at this frequency was presented for the first time in C13, 
and so it is not shown again here. The data were collected in a 
total bandwidth of 33.33 MHz split into 512 spectral channels. 
The flux density calibrator used was 3C 286, for which a flux 
density of 14.92 Jy was set using the Perley-Butler 2010 VLA 
valued]. Regular observations of the source 1822-096 were used 
for phase and bandpass calibration. Reduction and imaging were 
carried out with the NRAO Astronomical Image Processing Sys¬ 
tem (AIPS) package. The data from each day were fully reduced 
separately and then combined in a single data set. The reduction 
of the data consisted of the removal of radio frequency inter¬ 
ference (RFI) by hand rather than with the automatic flagging 
routines, and bandpass and gain calibrations. After initial cali¬ 
bration, the data were averaged in terms of frequency by col¬ 
lapsing the bandwidth to a number of 100 spectral channels. 
For the concatenated 624 MHz data set, we employed wide- 
field imaging based on a pseudo-three-dimensional multi-facet 
algorithm. In addition, we used a multi-scale CLEAN algorithm 
in AIPS, with four different scale sizes. Several self-calibration 
and imaging iterations were made to obtain the final image. 
The resulting GMRT image at 624 MHz has an angular reso¬ 
lution of 7".05 x 4".55, PA=44°.42, and an rms noise level of 
0.15 mJy beam -1 . We note the significantly high sensitivity (a 
factor of 10 better than the published image at 330 MHz) of this 
image at 624 MHz, being reported for the first time. 

On the other hand, the interstellar medium in the direction 
of G15.4+0.1 was investigated using data from the Southern 
Galactic Plane Survey in the 21 cm e miss ion line produced by 
HI (SGPS, iMcClure-Griffiths et al.1120051) . The SGPS images, 
combining data from the Australia Telescope Compact Array 
(ATCA) and the Parkes 64 m single-dish telescope, are valu¬ 
able in tracing large scale HI structures with an angular reso¬ 
lution of 3'.3 X 2'.1, PA=135°, and a velocity channel spacing 
of 0.82 km s -1 . 

2.2. Radio time-series observations 

Time-series observations towards the reported centroid of th e 
HESS source, / = 15°.41, b = 0°.17 dHofverberg et al.ll201lh . 
were carried out with the GMRT at 1404 MHz on May 11 and 
12, 2012, and at 624 MHz on June 17 and 18, 2012 (Project code 
22-015), simultaneously with the imaging observations. At each 
frequency, the GMRT was simultaneously used in two different 
modes to obtain two different time-series. In one mode, the out¬ 
put of 26 antennas were combined in an incoherent array (IA) 

1 For more detailed information regarding the flux density calibrators 
at VLA see the explanation of AIPS task SETJY or the website at 
http://www.nrao.edu/doc/vla/html/calib.shtml 


with an effective beam of 43' and 24' at 624 and 1404 MHz, re¬ 
spectively. The five X -ray point sources, identif ied in the XMM- 
Newton observations (lAbramowski et aljl2014lh all lie between 
2'-5' from the centre of this beam and were covered in the search 
for pulsation. A higher sensitivity phased array (PA) of 15 an¬ 
tennas, with an effective beam of 100" and 40" at 624 and 
1404 MHz, respectively, and pointed at the centroid of the HESS 
position, was used in the other mode. Both time-series were sam¬ 
pled every 61 ps with 512 channels across a 33.33 MHz band¬ 
pass at each frequency of observation. To accomplish the pe¬ 
riodic requirement of phase calibration for both imaging and 
phased array observations, the time-series were acquired in 13 
non-contiguous integrations of about 20 minute duration each at 
624 and 1404 MHz, respectively. In order to assess the sensitiv¬ 
ity of the search for radio pulsations, as well as the dispersion 
measures (DM) and the flux density estimates, we also carried 
out time-series observations in direction of three known pulsars: 
B1642-03, B1937+21,and J1901-0906. 

2.3. X-rays observations 

With the aim of performing a multiwavelength analysis of the 
emission from G15.4+0.1, we reprocessed data extracted from 
the on-line archive of the XMM-Newton observations^. The ob¬ 
servations in the direction of the remnant were made on October 
10, 2012, during revolution 2351 (Obs-ID: 0691390101). 

The data were reduced using the XMM-Newton Science 
Analysis System (SAS) software package Version 13.5.0. We 
applied standard processing to this data set to obtain clean event 
files. In order to mitigate the impact of high background flare 
activity, an appropriate screening was applied, extracting a light 
curve for photons above 10 keV, for the entire field of view. We 
selected events with FLAG = 0 and PATTERN <12 for both EPIC 
MOS1 and MOS2 cameras dTurner et al.ll200ll) . The resulting 
total effective exposure time of observation with the MOS1 and 
MOS2 cameras was 30 ks. 

In order to reveal the diffuse emission, point-like sources 
within the field of view were extracted by using the dmfilth 
tool of the CIAO (v. 4.6) reduction packageQ with POISSON 
statistics to fill the excluded regions. Finally, a single combined 
image in the energy range 1-8 keV from both cameras was ob¬ 
tained, which was then smoothed with a Gaussian kernel to a 
resolution of ~1'.4. 

3. Results 

3.1. Total intensity image of G15.4+0.1 at 624 MHz 

Figure [I] shows the new image of SNR G15.4+0.1 obtained at 
624 MHz using the GMRT. It constitutes the first ever high- 
resolution and sensitivity image of this source at low radio fre¬ 
quencies. The shortest baseline available in the GMRT observa¬ 
tions at 624 MHz is ~ 197/1, which implies that spatial scales of 
the G15.4+0.1 emission up to ~17.4 arcmin are well sampled in 
the resulting image. The radio emission from this remnant has 
the appearance of an irregular shell somewhat elongated in the 
Galactic north-east to south-wesQ direction (i.e. from upper left 
to lower right of Fig.[T} with an average diameter of 14.5 arcmin 

2 XMM-Newton Science Archive (XSA) 

http://xmm.esac.esa.int/xsa/index.shtml 

3 More information available at http://cxc.harvard, edu/ciao/ 

4 Although we do not utilize the equatorial coordinate system, we refer 
hereafter to north and south to set out locations around the Galactic 
plane. 


Article number, page 2 of 111 1 












L. Supan et al.: A complete radio study of SNR G15.4+0.1 from new GMRT observations 


(i.e. ~20 pc at the distance of 4.8 kpc derived in C13). The same 
overall total intensity morphology is also observed in the im¬ 
ages ofthisremnaiUa^30_andJ420^Hzje£ortedb^Brogan 
et al. d2006 ) and in C13 based on VLA and GMRT observations, 
respectively. At 624 MHz, the emission is brighter and patchy 
along the Galactic northwestern rim and from the northeastern to 
the central part of the radio shell, overlapping the region where 
shocked CO gas has been identified (C13); see also Fig. [4] in cur¬ 
rent work. The bright spot at l = 15°.52, b = 0°.19 corresponds 
to the 624 MHz radio counter part of the bipolar bubble HII re¬ 
gion G 0J^520+ 0.188 (Anderson et al.112011 1) , for which Supan 
et al. (l2014t) have determined a distance of 3.0 ± 0.3 kpc based 
on the new GMRT observations and HI emission line data. 

The radio emission becomes fainter and uniform towards the 
southern part of G15.4+0.1. Unlike the north edge, no sharp 
boundary is apparent in the south limb of the remnant. A gap 
in the radio emission is also evident along the western edge of 
the remnant. 



15.50 15.40 

Galactic Longitude (deg) 


15.30 


Fig. 1 . Total intensity image of SNR G15.5+0.1 at 624 MHz obtained 
with the GMRT. It does not include primary beam correction. The image 
was smoothed to an angular resolution of 10", resulting in an rms noise 
level of 0.48 mjy beam -1 . The unit of the values in the colour bar is 
mJy beam -1 . 


To accurately determine the flux density integrated over the 
whole 624 MHz emission from G15.4+0.1, we subtracted the 
contribution of the aforementioned HII region and the bright 
knot located at Z~15°.48, b~ 0°.20, which is likely to be an extra- 
galactic source, previously identified as NVSS J181806—152321 
(NED database^- The integrated flux density for G15.4+0.1 
based on the new GMRT observations at 624 MHz is 8.0 +1.1 Jy 
using the lPerlev & Butler! (i2013t) flux density scale. The error in 
our estimate reflects the uncertainty both in the determination of 
background emission and in the selection of the area over which 
we integrate the remnant flux. 


5 NED (the NASA/IPAC Extragalactic Database) is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under con¬ 
tract with the National Aeronautics and Space Administration. 


3.2. Pulsar search 

The time-series data for both the incoherent and coherent array 
for each frequency were analyzed using SIGPROC0 pulsar data 
analysis software using a 16 core 32-node High Performance 
Computing cluster (HPC) at the National Centre for Radio As¬ 
trophysics (NCRA, India). First, the channels affected by RFI 
were flagged from the data. The usable bandwidth of the obser¬ 
vations was somewhat less than 33.33 MHz owing to the pres¬ 
ence of narrowband RFI channels, which were removed, and was 
different for incoherent and phased array observations. These 
were then dedispersed to 1792 trial dispersion measure^ (DMs) 
covering a range of 0-815 pc cm -3 at 624 MHz. The DMs were 
spaced more coarsely at 1404 MHz, owing to a small amount 
of dispersion smear across each channel. Data were dedispersed 
to 360 trial DMs at this frequency covering 0-1200 pc cm -3 . 
A harmonic search was performed on the dedispersed time- 
series for each DM. Because of the non-contiguous nature of 
the time-series, an incoherent addition of spectra from individ¬ 
ual time-series segments was carried out before the harmonic 
search. Known RFI periodicities, such as power-line frequencies 
at 50 and 100 Hz, were excised from the resultant spectra before 
the harmonic search. The resultant candidate periodicities were 
carefully examined using a visual representation of the search 
results as shown in Fig. [2] The search typically produced 60 to 
70 candidates at each frequency for the phased array observa¬ 
tions, while the number of candidates for the incoherent array 
were larger. No pulsations above a threshold signal-to-noise ra¬ 
tio (S/N) of 8 were detected in either phased array or incoherent 
array observations at both frequencies. 

The sensitivity of the pulsation search was assessed us¬ 
ing observations of the three known pulsars, PSRs B1642-03, 
B1937+21, and J1901-0906, as mentioned in Section 12721 A 
procedure similar to that outlined in ICanellas et al . : d2012h was 
used for this purpose. The sensitivity of pulsation search depends 
on the duty cycle (ratio of pulse width to pulsar period) of the 
pulsar. As SNR G15.4+0.1 lies in the direction of the Galactic 
centre, any associated pulsar is likely to be scatter-broadened, 
particularly at 624 MHz, which affects the d uty cy cle of the 
pulsar. For example, using NE2001 models (iCordes & Laziol 
|2002|) with a distance of 4.8 kpc to the SNR (Cl3), we get an 
estimate of temporal broadening of about 13 ms at 624 MHz, 
whereas this value is about 300 /zs at 1404 MHz. Hence, we 
have plotted our sensitivity curves as a function of duty cycle 
in Fig. [3 to provide a more meaningful comparison in the pres¬ 
ence of temporal broadening. The individual curves correspond 
to the incoherent array and the phased array observations at 624 
and 1404 MHz, respectively. These curves provide 8 standard 
deviation upper limits, for pulsed emission with 10% duty cy¬ 
cle, of 250 and 300 /By in the phased array observations at 624 
and 1404 MHz, respectively, from any putative pulsar associated 
with the SNR G15.4+0.1/HESS J1818-154 system. Upper lim¬ 
its of 700 and 500 /By in the incoherent array observations at 
624 and 1404 MHz, respectively, are also implied for the five 
identified X-ray point sources within the extent of the SNR. We 
note that the incoherent array beam covers the entire remnant. 

3.3. A multiwavelength view of SNR G15.4+0.1 

We have explored the spatial co-existence of synchrotron ra¬ 
dio features and sites of emission at infrared and X-ray wave- 

6 http://sigproc.sourceforge.net 

7 The integrated electron column density along the line of sight in units 
of pc cm - ’. 
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Fig. 2. Visual representation of a pulsar search candidate. The plot shows the detection of the known pulsar PSR J1901-0906, following the 
procedure described in the text. The bottom left panel shows the peak in the S/N curve over a large range of dispersion measures (DM) values, 
whereas the top left panel indicates the peak in the S/N with a refined period and DM. The top right plot shows a greyscale representation of 
intensity as a function of pulse phase in abscissa and the bandpass in ordinate at the refined period and DM, while the bottom right plot shows the 
averaged profile of the pulsar. Two rotations are shown for clarity in this plot. This representation helps distinguish mistaken RFI candidates from 
genuine candidate pulsars. 



Fig. 3. Sensitivity curves for the pulsational search. The plot depicts the 
8 standard deviation upper limits obtained from our observations for in¬ 
coherent array at 1404 MHz (solid black), phased array at 1404 MHz 
(dashed-dot black), incoherent array at 624 MHz (solid blue), and 
phased array at 624 MHz (dashed-dot blue) as a function of duty cy¬ 
cle for a possible pulsar. 


lengths. In order to compare radio and far-infrared (FIR) mor¬ 
phologies, we made use of archival data at 70 and 250 pm ac¬ 
quired using the S PIRE i nstrument on Herschel Space Obsen’a- 
fonFI dGriffin et al.1120 Hi) . Far-infrared images are a very useful 
tool with which to delineate shock heated dust, particularly in 
regions where the SNR shock front is hitting dense interstellar 
medium. The spatial correspondence between the new radio im¬ 
age at 624 MHz and the FIR images is shown in Fig. [4}t; yel¬ 
low contours are included to help locate the molecular cloud 
discovered by Cl3. To perform the radio/IR comparison, the 
70 pm data were convolved to the same angular resolution as 
the 250 pm map (18")- From Fig. [4}i, the correlation between 
the CO molecular cloud and the spatial distribution of the emis- 

8 Herschel is an ESA space observatory with science instruments pro¬ 
vided by European-led Principal Investigator consortia and with impor¬ 
tant participation from NASA. 


sion at IR wavelengths is clearly evident. Indeed, the elongated 
structures that appear in the FIR images is completely bounded 
within the dense CO structure detected at velocities between 46 
and 50.3 km s 1 (C13). Moreover, significant enhancements in 
the dust emission at the two wavebands are located inside the two 
molecular clumps (referred to in C13 as clumps A and B). Such 
coupling between the dust emission in the FIR and the CO line 
emission provides additional support to the presence of dense 
molecular gas associated with the SNR. We thus interpret the in¬ 
frared emission as marking the SNR-molecular cloud interface 
in G15.4+0.1. Evidence of a similar situation has also been ob¬ 
served in the northern shock front of Cas A, in the region where 
the remnant encounters_a density enhancement in the interstellar 
medium (ISM) (iHines et al.ll2004l) . 

To obtain a more complete description of the infrared dis¬ 
tribution in G15.4+0.1, we also used data at 3.6, 4.5, 5.8, and 
8 pm taken with the IRAC camera aboard the Spitzer Space 
Telescope along with MIPSGAL images at 24 pm. Aside from 
the striking correspondence observed at FIR wavebands with the 
CO cloud, no discernible emission was detected in the near- or 
mid-infrared wavebands, which correlates with either the radio 
surface brightness or the proposed PWN in X-rays (the com¬ 
parisons are not shown here). This result contrasts with that ob¬ 
served in several PWNe in the Galaxy, for which there is a spatial 
agreement between synchrotron radio. X-rays, and near-infrared 
emission distribution (see e.g. the mult iwav elengt h vie w of the 
PWNe associated with SNRs 3C 58, ISIane etaf ( 20081) ; Crab 
ISandberg & Sollermanll 20091): G21,5-0 .9, Zaiczvket al J (120 1 2l) : 


and G292.0+1.8. iZharikov et all d2013l) ). 


The comparison o f the emission in radio an d X-ray bands 
was first presented by lAbramowski et al' .! t !2014l) using the im¬ 
age of the remnant taken from the 330 MHz VLA survey of the 
Galactic plane (GPS. iBrogan et ak l [2005al ) and XMM-Newton 
observations. Here, we redo this analysis in the light of the new 
GMRT data by combining the synchrotron emission at 624 MHz 
from the source, which shows a highly structured interior com¬ 
pared to that observed at the lower radio frequency, with XMM- 
Newton reprocessed data obtained from X-ray energies between 
1.0 and 8.0 keV. The radio/X-ray morphological correlation is 
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shown in Fig. IS>. Contours for the molecular gas interacting 
with the SNR (Cl3) together with the region where y-ray emis¬ 
sion was observed are also presented in the radio/X-ray compar¬ 
ison. Fig.UJt reveals little correspondence between both emitting 
plasmas. The X-rays associated with G15.4+0.1 do not show 
evidence of a limb-brightened X-ray shell, but most of the X- 
ray emission appears broadly contained within the radio rims 
of G15.4+0.1. Moreover, the comparison shows X-ray emission 
extending slightly beyond the projected radio boundaries of the 
remnant. In addition, weak and diffuse X-ray emission is present 
filling the spatial gap observed in radio all along the western limb 
of the remnant. 

Figure [4}> also demonstrates a significant lack of correspon¬ 
dence between the X-rays and the molecular material, delineated 
by yellow contours. This result may not be surprising if obscura¬ 
tion due to the molecular cloud detected between the observer 
and the remnant occurs. On the contrary, the morphology in 
X-rays overlaps part of the y-ray emission area delineated by 
the white circle in Fig. [4b . Such a X-ray/y ray correlation led 
lAbramowski et al] (120 14t) to the conclusion that the X-ray emis¬ 
sion is the counterpart to the TeV source originating in a pulsar 
wind nebula of a yet undiscovered pulsar. This interpretation, 
however, may not have a robust observational support since the 
authors remark that both thermal and non-thermal models are 
equally satisfactory to fit the spectrum of the diffuse X-rays in 
G15.4+0.1, which made somewhat uncertain the true origin of 
the detected X-ray emission. Finally, we found no radio coun¬ 
terparts corresponding to the five X-rays point sources detected 
using XM M-Newton observations a t the location of the emission 
in y rays dAbramowski et al.ll2014t) . 

4. Radio spectral analysis in SNR G15.4+0.1 

4.1. The integrated spectrum 

To determine the global spectral index (S oc y") of G15.4+0.1, 
we have used the new flux density estimate at 624 MHz together 
with that obtained from GMRT data at 1420 MHz (Cl3) and 
measurements at 330, 1400, and 2700 MHz as determined from 
obse rvati ons carried out with the VLA and Effelsberg telescopes 
iBrogan et al.ll200 6). In Table Q] we summarized the flux densi¬ 
ties integrated over the whole extension of this remnant at the 
mentioned frequencies. The listed values constitute all the mea¬ 
surements published up to present for the source. A correction 
factor was applied to set all the values to the iPerlev & Butleil 
(12013h absolute flux density scale. In particular, the correction to 
the flux density estimates at 330 and 624 MHz was made using 
the low frequency extension proposed by those authors. Figure [5] 
shows the resulting radio continuum spectrum of the remnant; 
the determinations from GMRT observations are plotted as open 
pentagons symbols. A weighted fit to all the available flux densi¬ 
ties yields a spectral index a=- 0.62 ± 0.03, which i s comp atible 
with the value derived previously by [Brogan et al.l (2006). This 
spectral behaviour is co nsistent with that measured in other shell 
type SNRs dKothes et al.1120061) . 

4.2. Spatial variations in the radio spectral index 

We have investigated the spectral changes as a function of the 
position within the remnant by using the image at 624 MHz and 
the one at 330 MHz (see panel (f) in Fig. [6} kindly provided by 
C. Brogan (iBrogan et al.ll2006b . The analysis was done by con- 
stru cting tomographic mans between them ( Katz-Stone & Rud- 
nick ll997h . " 


Table 1. Integrated flux densities for G15.4+0.1. 


Frequency 

(MHz) 

Scaled 

flux density (Jy) 

References 

330 

11.3 + 0.3 

Brogan et al. (2006) 

624 

8.0+ 1.1 

This work 

1400 

4.7 ±0.8 

Brogan et al. (2006) 

1420 

4.7 ±0.2 

C13 

2700 

2.9 ±0.3 

Brogan et al. (2006) 



Log t/(MHz) 

Fig. 5. Integrated radio continuum spectrum of SNR G15.4+0.1. The 
open symbols correspond to GMRT flux densities measureme nts at 624 
and 1404 MHz. The rest of the data were taken from IBrogan et all 
120061k All the_data_were i ^ied -i to -- a ii common -i flux_densit^^cale_(Per- 
ley & Butler l2013f) . The solid line shows the best fit obtained with a 
single power law to all the values, which yields a global spectral index 
a=-0.62 ± 0.03. 

Basically, the procedure involves the generation of a se¬ 
ries of residual images S, according to the expression S t = 
S 330 - (pf) 1 S 624 , where S 330 and S 624 are the maps at 330 and 
624 MHz, respectively, and a t represents a trial spectral index 
value. If the assumed spectral index coincides with the “true” 
value, one should expect the S t image to be equal to zero. On 
the contrary, features in the residual map that appear positive 
(negative) will be associated with local spectral indices that are 
steeper (flatter) than the assumed cr t value. The technique has the 
advantage of delineating spectral features that can overlap along 
the line of sight while avoiding any zero level dependence. 

In Fig. [6] we reproduce the gallery of tomographic maps be¬ 
tween 330 and 624 MHz. To construct them, the radio images 
were clipped at the 3cr level of their respective noise levels. In 
order to obtain a reliable determination of the spectral proper¬ 
ties, the range of the spatial scales measured at 624 MHz was 
matched in the t/v-coverage to that of the 330 MHz data. Ad¬ 
ditionally, the images were aligned and interpolated to identical 
projections to avoid positional offsets. The a t value shown in 
the maps varies from +0.1 to -0.7. For the sake of illustration 
the 330 MHz total intensity image of G15.4+0.1 is shown in the 
bottom right frame in Fig. [6] 

Spatial variations in the spectral index are recognizable 
across the synchrotron emission from G15.4+0.1. The central 
part of the remnant shows components of spectral index flat¬ 
ter (a varying between —(-0.1 and -0.5) than the values in the 
ridge. We note that the very flat spectral component observed as 
a dark area in the southeastern border of the remnant (l~ 15°.40, 
b~ 0°.15), within which the spectral index appears to have posi- 
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Fig. 4. Multiwavelength view of SNR G15.4+0.1. a) A color composite image showing the new radio image at 624 MHz of G15.4+0.1 (in blue) 
and FIR emission as observed by Herschel at 70 pm (in red), and 250 /jm (in green); b) Comparison of the radio continuum emission at 624 MHz 
of the SNR (in blue) and the X-ray emission (in yellow) as observed by XMM-Newton in the 1-8 keV energy range. The dashed white circle 
indicates the size and location of the TeV source HESS J1818— 154. In both panels, the yellow contours delineate the 13 CO line emission in the 
46-50.3 km s _1 range (C13). 


tive values (see panel corresponding to a,-+(). 1 in Fig. [6]>, cor¬ 
responds to a bordering clipped area in the map and should not 
be considered a real spectral component owing to the decrease 
in the signal-to-noise ratio observed in this region in the VLA 
330 MHz image in comparison with the image at 624 MHz. In 
general, the central flattening traced in the 330/624 MHz com¬ 
parison matches the zone where the interaction with the dis¬ 
covered CO molecular cloud is taking place (see Fig. [4]». Vari¬ 
ations in the spectrum are generally expected in regions where 
the SN shock is impacting high density interstellar material as a 
consequence of the density increase and the deceleration of the 
shock in the interacting regions. For instance, spectral flattening 
were found between 74 and 330 MHz for the two well-known 
SNRs/molecular clouds systems marking the ionized interface 
of molecular clouds interacting w ith the shock front in 3C 391 
dBrogan et al.ll2005bl) and IC 443 (iCastelletti et al.ll201 ll) . 

The tomographies also show a conspicuous bullet-like fea¬ 
ture standing out as a bright spectral component (i.e. positive in 
the St maps) located at the northeast part of the shell (/~15°.48, 
b~ 0°.25). It is characterized by a spectral index even steeper than 
-0.7 because it appears as a positive residual in all tomographic 
images. We found no catalogued extragalactic objects at the po¬ 
sition of this feature. 

Along the westernmost edge (left side) of the SNR, there is 
an indication of a positive correlation between the local spec¬ 
tral index and the total radio intensity features in the sense that 
the bright ridge tends to be associated with steeper indices. A 
similar situation was found in young SN Rs (see e.g. the analy¬ 
sis^ SNRs G39.2-0.3 and G41.1-0.3 in I Anderson & Rudnicki 
119931) . In those cases, it was suggested that the radio brightness 
and spectra are being regulated by different mechanisms. In addi¬ 
tion to this behaviour, most of the structures in the tomographies 
traced at (+=-0.5 and a t =-0.7 disappear against the background 
indicating that the spectral behaviour can be characterized by a 
mean spectral index —0.6, which is consistent with the inte¬ 
grated spectral index of this SNR derived by fitting the total flux 
densities from 330 to 2700 MHz (see Sect. 4.1). Overall, the ob¬ 
served radio continuum spectral distribution across G15.4+0.1 is 
compatible with what is expected in the framework of diffusive 
shock acceleration models. 


5. Neutral gas distribution around G15.4+0.1 

5.1. HI brightness morphology and physical parameters 

We explored the interstellar medium around G15.4+0.1 in the 
21 cm line to distinguish some morphological correlations be¬ 
tween the distribution of the neutral gas and the remnant that 
might account for the SNR shock effects on the surrounding 
medium. 

Figure [7] includes velocity channels maps from the SGPS 
data cube covering an interval from 41.6 to ~ 108 km s _1 . All 
the velocities are measured with respect to the local standard of 
rest (LSR). Each frame was obtained by integrating the data cube 
every 8.2 km s . An appropriate background level equal to the 
mean value of each integrated image was subtracted from every 
map. A contour line delineating the outer boundary of the radio 
continuum emission at 624 MHz from G15.4+0.1 was overlaid 
to the HI distribution for ease of comparison. 

The distribution of the neutral gas in the velocity range 
from 41.6 to 50 km s 1 only shows bright clouds spread over 
the field without any correspondence with the radio continuum 
shell of G15.4+0.1. A similar behaviour is found between 0 and 
41.6 km s and therefore, the images are not shown. In the im¬ 
ages corresponding to velocities between 50 and 91.1 km s , 
we observed the best morphological signature of neutral gas re¬ 
lated to the SNR. Indeed, at ~54 km s _1 , the HI distribution ap¬ 
pears to be distorted at the position of the synchrotron emission 
from G15.4+0.1. Moreover, it is in the 58.1-91.1 km s range 
where we see a cavity in the HI emission that closely matches 
G15.4+0.1. In particular, at velocities from 58.1 to 74.6 km s , 
the HI distribution has an almost complete shell-like morphology 
along the eastern and western peripheries (left and right sides) of 
the remnant, which looks broken toward the northeastern border. 
The presence of this shell is compatible with the central sys¬ 
temic velocity v sys = 60 km s ' derived in C13 for the SNR. The 
expanding shock wave of the SNR on the HI surroundings might 
be responsible for the creation of the observed HI structure. An 
alternative origin for the HI shell may be a pre-existing stellar 
wind bubble formed for instance by a group of OB stars, one of 
which could have been the SN precursor of G15.4+0.1. In this 
particular scenario, the brightest synchrotron emission observed 
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Fig. 6. Tomographic images for SNR G15.4+0.1 constructed from data at 330 and 624 MHz. The data at both frequencies were matched to the 
same uv-coverage. The radio spectral index a t is indicated at the top right corner of the panels. The grey colour bar in the middle column represents 
the scale of the residuals components, which is the same for all tomographies. In the greyscale representation, spectral components for which a is 
steeper than the assumed a t are shown as light regions (positive residuals), while flatter components are seen as dark features (negative residuals). 
A contour delineating the radio continuum emission at 330 MHz is included for reference. The total intensity image of G15.4+0.1 at 330 MHz is 
plotted in the bottom right frame ( IBrogan et al.ll2006li . 


towards the northwest face of the remnant could be where the 
shock is interacting with the inner wall of the bubble. Further¬ 
more, the fact that 13 CO emission was detected on the shell’s 
central face of G15.4+0.1 at approximately 50 km s~ ] (C13) 
may indicate that the molecular cloud is outside the front wall 
of atomic hydrogen; perhaps the inner wall is also the place 
where the synchrotron emission originates. Examples of bub¬ 
bles created by the action of the precursor’s stellar winds, inside 
of which evolves a SNR shock wave, are those discovered sur¬ 
rounding the remnants Kes 79 dGiacani_e^ - al. 2009), 3C 434.1 
dFoster et alJl2004 . G106.3+2.7 dKothes et alJl200lh . In any 
case, whether the HI shell is formed by neutral gas accelerated 
by the SNR shock or by stellar winds from the SN progenitor, the 
X-ray emission detected inside the radio shell of the SNR rein¬ 
forces the idea that the center has been evacuated. It is also worth 
noting that at higher velocities there is no evidence of the natu¬ 
rally expected cap-like features associated with an HI shell. The 
non-detection of possible caps does not imply that these struc¬ 
tures do not exist. They could be undetected because of confu¬ 
sion caused by unrelated foreground and background inhomo¬ 
geneities in the ISM. The caps, produced by accelerated HI gas, 
have been observed either in e mission as central conc entrations 
projected within the remnant (IVelazquez et al.l 12002 1) or in HI 
sclf- absorntion in an expanding stellar wind bubble wall ( Foster 
et al. l2004l) ~ At velocities higher than —91.1 km s -1 , the 21 cm 
line emission does not correlate with the SNR, and neutral gas is 
preferentially oriented parallel to the Galactic plane. 


In what follows, we derive a number of parameters charac¬ 
terizing the HI shell. For this structure, centered at l~ 15°.42, 
b~ 0°.16, we assigned a radius of ~9'.8 or —13.7 pc at a dis¬ 
tance of 4.8 kpc. If we assume that the HI emission is optically 
thin, the integration for velocities between 50 and 91 km s , 
for which we observed the most evident signatures of associa¬ 


tion between the HI gas and the SNR, yields an average column 
density Ah — 4 x 10 20 cm~ 2 . The total mass of atomic gas that 
forms the shell is about 1.9 x 10 3 M e . This value is in fact an 
upper limit due to the impossibility of disentangling related and 
unrelated HI gas. If the HI shell was formed by ISM swept up by 
the SNR shock and assuming that before the SN event the total 
neutral mass was uniformly distributed inside the volume of the 
HI shell, the derived mean density of hydrogen nuclei in the am¬ 
bient environment is rjo ~7 cm -3 , which is larger than the typical 
interstellar hydrogen densities of ~1 cm~ 3 (averaged over the 
cold, warm Jii andho^ i gas i £hase_constituents_ofthe i JSM fc McKee 
& Ostriker Tl977l) . Our result corresponds to a reasonable kinetic 
energy supplied by the SN explosion into the the surrounding 
medium of about 7.6 x 10 48 erg. This value was estimated by 
adopting an expansion velocity for the HI shell of ~20 km s _1 , 
based on the systemic velocity of the neutral gas derived in Cl3 
(v sys ~60 km s 1 ) and the Hl-continuum correspondences seen 
between ~50 and 91 km s -1 . Taking into account errors in the 
selection of the boundaries for integration, confusion from back¬ 
ground or foreground unrelated HI gas, as well as uncertainties 
in the determination of the distance, we derived the parameters 
for the HI shell with a mean relative error of 45%. 

As a by-product, on the basis of the new radio continuum im¬ 
age at 624 MHz and the results derived from the HI data, we cal¬ 
culated physical parameters of the remnant G15.4+0.1. By com¬ 
bining the ambient density 770 ~1 cm~ 3 , the radius of G15.4+0.1 
Rsnr~7'.2 (—10.1 pc at the distance of 4.8 kpc) with the velocity 
of the shock, v S h, we can roughly calculate the energy released 
by the SN explosion to the ISM from the expression given in the 
model discussed bv IChevalieil (1 974 ). 


(1) 
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To estimate the shock velocity, we used v s h = 
66.5 (§ff) -2 ' 23 km s _1 and obtained v S h = 370 km s _1 for 
a shock radius set to /?snr —10.1 pc. We understand that this 
value represents a reliable bound for the velocity of the SN 
shock before impacting the dense molecular cloud located in 
front of G15.4+0.1; the velocity range of 46-50.3 km s _1 for 
the cloud corresponds to a distance of ~4.2 kpc (C13). To test 
the robustness of our interpretation, we also estimated the shock 
velocity using the physical parameters derived for the 11 CO 
cloud interacting with SNR G15.4+0.1, as calculated with the 
expression v s h = v c i (^ci/t/int) 1 ^ 2 , where rj c \ and v c i are the density 
and the expansion velocity of cloud, respectively, and ;/i llt is the 
intercloud density assumed to be 1 cm 3 dDubner et al .11 19991) . 
In Cl3, it was demonstrated that the CO cloud interacting 
with G15.4+0.1 has a density of ~1.5 xlO 3 cm~ 3 and it is 
expanding at ~10 km s _1 . In these circumstances, we derived 
a shock velocity of v S h - 380 km s _1 , which agrees within the 
uncertainties of the method with that obtained using HI data. 
Finally, for an average shock velocity of 375 km s 1 and using 
the ambient density along with the radius of the shell derived 
from our study, we obtain an initial energy of ~2.6 xlO 51 erg, 
which agrees with the standard value. 

To account for the age of G15.4+0.1, we employed the 
dynamical evolution model of IChevalierl (1 1 974l) and obtained 
t a (§p|f) 3 ' 23 10 5 - 8.2 xlO 3 yr, where we used a radius for 
the SNR of Rsnr~ 10.1 pc at a dista nce of about 4.8 kpc. This 
result differs from that calculated by lAbramowski et all (l2014lf 
(2.5xl0 3 yr) based on a Sedov-Taylor and the assumption of an 
ambient density around G15.4+0.1 of 1 cm -3 . 

5.2. Modelling the broadband spectrum of G15.4+0.1 

The analysis in lAbramowski et all (120141) to determine the ori¬ 
gin of the y-ray emission in the region of G15.4+0.1 excludes a 
hadronic contribution. In modelling the spectral energy distribu¬ 
tion (SED), these_authors_used radio data at 330 and 1400 MHz 
from the GPS dBrogan et al.ll2005ah and the Mu lti-Array Galac¬ 
tic Plane Imaging Survey (MAGPIS. iHelfand et al.ll2006h . re¬ 
spectively, whose fluxes were treated as upper limits, X-ray 
data from XMM-Newton, and y-ray observations obtained with 
H.E.S.S. In this section, we first demonstrate that, with the avail¬ 
able data, a hadronic model in which the TeV y-ray flux is pro¬ 
duced through the interaction of accelerated protons with ambi¬ 
ent protons, mainly from the molecular cloud located foreground 
G15.4+0.1, is still compatible with the broadband spectrum of 
the G15.4+0.1/HESS J1818—154 system. In addition, in the light 
of the ne w GMRT data, we revisit ed the leptonic model pre¬ 
sented in lAbramowski et all (i2014f) to fit the multiwavelength 
spectral data. 

The energy spectrum of the accelerated particles (electrons 
and protons) in the emission region is described by a power law 
with an exponential cut-off. 


^ = K ( , P exp (~E/E cutep ) , (2) 

where the subscripts e and p refer to the particle species (elec¬ 
trons or protons), and E, K, T, and £ cut are the particle energy, 
a normalization constant, the spectral index, and the cut-off en¬ 
ergy, respectively. We refer to lAharonian et al. (2010 ) and refer¬ 
ences therein for calculation details of the synchrotron emissiv- 
ity produced by electrons. Regarding the inverse Compton (IC) 
and the non-thermal Bremsstrahlung emission, we followed the 


method given in Jones (11968 ) for the former, while the method 
presented in lBaring et al.l (1 1999h was used to model the electron- 
elec tron Bremsstrahlung i nteraction and in lKoch & Motzl (119591) 
and ISturner et al.l (Il997h for the electron-ion Bremsstrahlung 
process. On the other hand, we determine the y-ray spectrum 
due to the neutral gion de cays b y using the ppfrag program 
(IKachelrie s & Qstanchenkoll2012 ). In the calculations, the low 
and high energy parts of the spectrum were separately consid¬ 
ered. For energies of the y rays smal ler than 50 GeV, we used 
the parametrization presented in lKamae et al . (1200 7). while for 
larger values of energy we based our calculations in the interac¬ 
tion model QGSJET-II-(kfl (IOstat>chenkol20lH) . 

The XMM-Newton observations show faint X-ray emission, 
about 9 arcmin in size, which partially coincides with the size 
of the HESS J1818—154 source (see Fig. Q. The origin of the 
X-rays was considered non-thermal and based on their spatial 
correspondence with the y rays, an asso ciation between them 
was suggested (lAbramowski et ahll2014l) . We model the spec¬ 
trum of G15.4+0.1 from radio to the TeV y rays by considering 
a non-thermal nature of the X-ray emission. This implies that 
the whole X-ray flux and part of the radio flux is produced by 
synchrotron radiation of electrons in the turbulent magnetic field 
in the region interior to t he SNR shell, the PWN proposed in 
lAbramowski et all (120141) . We refer to this zone as Region I. 
The remaining part of the radio flux is also generated by syn¬ 
chrotron radiation, but corresponding to a different population 
of electrons that are placed in the shell and referred to hereafter 
as Region II. Figure [8] shows the broadband spectrum for the 
SNR G15.4+0.1/HESS J1818—154 system. The radio data for 
Regi ons 1 and II include the observations presented in Brogan 
et al. d2006() and the new GMRT data at 624 and 1420 MHz. 
The spectrum also inclu des the best fit to th e X-ray data using 
a power law taken from lAbramowski et al~ J 20141), and upper 
limit y-ray fluxes obtained by Fermi -LAT ( Acero et al.ll2013h 
and H.E.S.S (lAbramowski et al.l2014l) . Although upper limits on 
the fluxes at GeV energies are represented, they were not used 
in fitting the spectrum as they are too far from HESS data. The 
best fit to the radio data corresponding to the interior SNR re¬ 
gion is shown by a dashed curve labelled I. The intensity of the 
magnetic field and the cut-off energy are fixed during the fit¬ 
ting procedure in order to be consistent with the power-law fit 
to the X-ray data. The normalization constant K and the spec¬ 
tral index T in Eq. (0 are taken as free fit parameters. The in¬ 
tensity of the magnetic field used is B l = 25 pG and the cut¬ 
off energy is log(£/!„, f /eV) = 12.83. These values are chosen 
in such a way that the inverse Compton component has a sub¬ 
dominant role at very high energies. It is worth mentioning that 
scenarios with larger values of the magnetic field intensity and 
smaller cut-off energies of the electron component, such that 
{Elute) 2 x B l = 1.1536 x 10 21 eV 2 G, also describe properly 
the X-ray data, making the inverse Compton component even 
smaller. The best-fit model for Region I yields a spectral index 
rj = 2.35 + 0.14. 


We also considered a contribution to the leptonic mechanism 
from accelerated electrons in the SNR shell. The radio flux, cor¬ 
responding to the shell region, is obtained by fitting the total 
radio flux, by using the best fit to the radio data corresponding 
to the internal region. In this case, we fixed the cut-off energy 
and the intensity of the magnetic field to log(£ c 13 , e /eV) = 11.6 
and B 11 - 62.5 pG, respectively. This choice of the parameters 
prevents the synchrotron flux that originated in the shell region 


9 The last version of the QGSJET model was updated with recent 
Large Hadron Collider data. 
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Fig. 7. Images of 21 cm line emission in the direction of SNR G15.4+0.1. Each map corresponds to the integration of the neutral gas every 
8.2 km s _1 . The range of velocities is indicated at the top right corner of the panels. The colour display was kept for all images between 633 
and 1210 K km s _1 . The 0.7 mjy beam -1 contour from the 35" resolution 624 MHz image is included in each panel to facilitate the comparison 
between the radio continuum emission from the remnant and its surroundings. 


from contributing to the X-ray part of the SED. We note that 
this combination of the two parameters is not unique. Smaller 
values of the cut-off energy, which are consistent with the data, 
produce smaller values of the high energy flux originated by the 
inverse Compton process. The spectral index obtained from the 
fit is F 31 = 2.16 + 0.08. It is evident that in our model the inverse 
Compton process alone cannot explain the emission at TeV en¬ 
ergies from Region I or II. Figure [8] also includes the spectrum 
corresponding to the Bremsstrahlung radiation. In this case, we 
assumed that only Region II has enough density to produce a 
non-negligible flux. The calculation was done by fixing the pro¬ 
ton density, both molecular and atomic, to the value n = 6.4 cm -3 
averaged over Region II. We note that, as in the IC model, it is 
possible to reduce even more the Bremsstrahlung effects dimin¬ 
ishing the value of the cut-off energy of the accelerated electrons 
of the shell region. 

It is worth mentioning that the magnetic field intensity of the 
two regions required to describe the data is much larger than the 
field corresponding to the interstella r me dium, for which its typi¬ 
cal strength is 1-2 /jG ( Foster et al.l2013l . and references therein). 
Thi s is typi cal in hadronic models of the y-ray emission (see e.g. 

ICardillo et al1l20l4 . 

We determined the hadronic contribution of the y-ray flux 
resulting from the neutral pion decay of accelerated protons (il¬ 


lustrated by a solid line labelled as n° in Fig. [8}. Since neither 
the H.E.S.S. data nor the Fermi-LAI upper limits are enough to 
constrain the spectral index of the proton component, we used 
the spectral index of the electron component fitted for the syn¬ 
chrotron emission in Region II, that is T 31 =2.16. Taking into ac¬ 
count the normalization constant and the cut-off energy as free 
fit parameters, we obtained for the proton component a cut-off 
energy log (E cu , tP /eV) = 14.6 + 0.5. 

On th e other hand, we revisited the one-zone model pre¬ 
sented in lAbramowski et afl (12014 1 replacing the upper lim¬ 
its on the radio fluxes considered in that work by reliable flux 
density measurements of the radio synchrotron emission in the 
interior of G15.4+0.1 (called Region I in the current work). 
This region is spatially coincident with the X-ray emission at¬ 
tributed by those authors to the same electron population re¬ 
sponsible for the y-ray emission through the IC process. In this 
case as well, a power law with an exponential cut-off was as¬ 
sumed to model the spectral energy distribution of electrons. As 
was done for the hadronic scenario, we imposed the condition 
(E cu t,e) 2 x B = 1.1536 X 10 21 eV 2 G, which ensures an ade¬ 
quate description of the X-ray data. We would like to note that 
also in this case the Fermi-LAI upper limits do not constrain the 
model. The corresponding spectrum is shown in Fig. [9] For this 
simple model, we obtain a spectral index F e = 2.41 ± 0.09, a 
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Fig. 8. Spectral energy distribution of the SNR G15.4+0.1. The filled 
squares and triangles with error bars correspond to the emission from 
the whole remnant and fr om th e interior of the shell (Region I), respec¬ 
tively (data from lBrogan et alfeood . and this work). The shadowed re¬ 
gion encompassing a solid line corresponds to the fitted power law and 
the one sigma region of th e X-ray data. The HESS T eY data points are 
indicated by filled circles! Abramowski et al. 20 14j), while the arrows 
poin ting do wnward s show the upper limits in the Fermi -LAT spectrum 
(iAcero et al.ll2013l) . The dashed lines labelled I and II mark the spectra 
from Region I, assuming non-thermal X-ray and radio emission origi¬ 
nating in the interior of G15.4+0.1, and from Region II in the shell of the 
remnant, respectively. The modelled spectra from the synchrotron radia¬ 
tion by electrons, inverse Compton on the CMB, Bremsstrahlung, and y 
rays produced by neutral pions generated in proton-proton interactions 
are indicated by the curves labelled S, IC, B, and ti 0 , respectively. 

magnetic field B = 3.9 ± 1.4 pG, and a cut-off energy of elec¬ 
trons log {E cute /cV) = 13.303 ± 0.071. The fit corresponding to 
this sc enario is as good as the one obtained in lAbramowski et alJ 
(120141) . However, the fitted parameters derived in this work are 
restricted to smaller intervals. This is due to the radio data points 
of the internal region used to fit the SED. 



Fig. 9. Broadband fit to radio (filled triangles) and y-ray observations 
obtained with H.E.S.S. (filled circles). Upper limits to the flux densities 
from Fenni-LAfT (arrows pointing downwards) and the best fit to the X- 
ray data (green hatching) are also included. The radio data correspond 
to the internal region of SNR G15.4+0.1. The curves labelled S and 
IC show the model of the synchrotron radiation and inverse Compton 
mechanism on the CMB by accelerated electron, respectively. 

The observational data is compatible with both scenarios an¬ 
alyzed here; more data is needed, especially in the energy region 
relevant to the Fermi-LAT observatory, in order to make progress 


in the understanding of the most important physical processes 
taking place in this source. In any case, it is important to empha¬ 
size that a hadronic scenario is compatible with the observational 
data. 

6. Conclusions 

We presented new full-synthesis imaging of the SNR G15.4+0.1 
obtained from observations at 624 MHz carried out with the 
GMRT. We measured an integrated flux density S - 8.0 +1.1 Jy 
at 624 MHz. Based on the combination of our estimate with 
those previously published, after bringing all values to the same 
absolute flux density scale, we derive a spectral index a——0.62 ± 
0.03 for the whole remnant. The estimated global spectral in- 
d ex is c onsistent with measurements for other typical SNR shells 
dKothes et al.ll2006i) . 

The new GMRT data were used to analyze the correspon¬ 
dences with FIR and X-ray observations towards the remnant. 
We found an impressive correlation between the infrared emis¬ 
sion detected by Herschel at 70 and 250 pm and the cloud of 
13 CO colliding with the SNR shock reported in C13. In other 
regions of G15.4+0.1, however, no IR counterpart is observed 
to the radio emission. On the other hand, from the comparison 
of the new 624 MHz image with re-processed XMM-Newton 
observations between 1 and 8 keV, we can confirm that the 
synchrotron radio emission surrounds the low surface bright¬ 
ness X-ray radiation. No counterparts were observed in radio to 
the five point-like sources detected in the XMM-Newton field. 
We also conducted a search for a neutron star associated with 
G15.4+0.1 through time-series observations performed with the 
GMRT, which provided negative results. 

The analysis of the spatial spectral index variations made be¬ 
tween existing observations of G15.4+0.1 at 330 and the new 
GMRT data at 624 MHz revealed a radio spectrum steepening 
from the weak interior to the brighter periphery of G15.4+0.1. 
A similar spectral behaviour is also found in other young SNRs. 
The overall distribution of the local spectral radio index could 
be described by diffusive shock acceleration by considering the 
effect of the ISM gas. We did not recognize any distinct fea¬ 
ture in the radio brightness at 624 MHz or in the 330/624 spec¬ 
tral distribution that could be considered a synchrotron nebula, 
powered by the wind of a yet undetected pulsar. Furthermore, 
our search for radio pulsations inside the synchrotron shell did 
not reveal a pulsar, with the upper limits on the mean flux den¬ 
sity of any associated putative pulsar of 250 and 300 pSy at 624 
and 1404 MHz, respectively, towards the centroid of the TeV 
source. In addition, the flux density upper limits that we esti¬ 
mated for any putative pulsar morphologically coincident with 
SNR G15. 4+0.1, including the five po int-like X-ray sources re¬ 
ported by lAbramowski et aD ( 2014 ). were of about 700 and 
500 /./Jy at 624 and 1404 MHz, respectively . 

Our study of the neutral ambient gas in the direction of 
G15.4+0.1 unveiled a large (about 13.7 pc in radius at a distance 
of 4.8 kpc) and incomplete HI shell-type structure. We deter¬ 
mined a mass for this shell of 1.9 x 10 3 M Q and a kinetic energy 
of 7.6 x 10 48 erg. On the basis of this study and the new SNR im¬ 
age at 624 MHz, we estimated that G15.4+0.1 was created as the 
result of 2.6 x 10 51 erg released in a SN explosion that occurred 
about 8.2 x 10 3 yr ago. 

Finally, we discussed different models to fit the broadband 
emission from the SNR G15.4+0.1/HESS J1818-154 system. 
In particular, we demonstrated that taking into consideration the 
multiwavelength data known to date, a purely hadronic picture 
is sufficient to account for the y-ray data. This is in contrast to 
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the explanation given in lAbramowski et ahi (l2014t) for which the 
X-ray and TeV emission have the same origin. We recognize 
that more comprehensive X-ray and y-ray data (especially in the 
MeV-GeV band) are required to improve the constraints on the 
observed SED of the system and distinguish between leptonic 
and hadronic explanations for the origin of HESS J1818—154. 
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